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Abstract:  This  work  introduced  the  theory  and  simulation 
results  of  employing  the  Kalman  filter  in  an  airborne  fire 
control  radar  tracking  system.  It  discussed  problems  such  as 
target  dynamics  modeling,  linear  filtering  and  simulation,  linear 
filtering  approximation,  sensitivity  simulation,  maneuvering 
target  tracking  and  adaptivity.  Simulation  results  showed  that 
the  system  performed  far  better  than  typical  systems.  Further¬ 
more,  its  adaptivity  to  maneuvering  targets  was  better. 

Abstract:  An  application  of-  tha  Kalman  filtering  theory  and  simulation  results  with  com¬ 
puter  are  developed  in  tracking  system  of  airborne  firecontrol  radar.  The  dynamic  mathe¬ 
matical  model  of  target,  the  linear  filtering  and  simulation,  approximation  by  the  linear  filter¬ 
ing.  sensitivity  simulation,  maneuvering  target  tracking  and  adaptive  filtering  and  etc  are 
showed.  Simulating  results  show  that  this  system  is  better  than  general  system  on  the  perfor¬ 
mance  and  the  adaptive  capability  of  maneuvering  targets  greatly. 

I.  Introduction 

Since  the  introduction  of  the  Kalman  filter,  it  has 
primarily  been  used  in  tracking  and  processing  of  targets  in 
orbits  and  suborbits  or  during  re-entry.  However,  very  little 
research  on  piloted  moving  targets  has  been  conducted.  Since  the 
seventies,  a  great  deal  of  effort  was  devoted  to  the  tracking  of 
maneuvering  targets.  The  results,  however,  are  usually  limited 
to  specific  targets,  detectors  and  environments.  Hence,  there  is 
a  need  to  find  a  simple  and  appropriate  model  to  approach  a 
moving  target.  The  Kalman  filter  was  used  to  build  an  advanced 
tracking  system.  This  work  absorbed  the  valuable  portions  in 
References  [1-4]  and  obtained  specific  results  through 

♦completed  in  March  1982,  revised  in  March  1983. 
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theoretical  analysis  and  simulation.  It  provided  a  basis  for 
engineering  design. 

The  single  pulse  (or  conical  scanning)  radar  tracking  system 

discussed  in  this  work  is  shown  in  Figure  1.  The  pitch  angle 

error  signal  Eg,  azimuthal  angle  error  signal  E^,  range  error 

signal  Er  and  speed  error  signal  Ev  were  expressed  digitally  by 

the  radar  receiver  and  signal  processor.  The  sampling  frequency 

of  these  data  could  reach  several  dozen  kHz.  The  recurrence 

cycle  of  a  Kalman,  however,  is  usually  10^  Hz  or  slightly  higher 

to  satisfy  the  requirements  in  closed  loop  tracking.  Therefore, 

in  order  to  minimize  loss  of  information,  a  prebalancer  was  used. 

To  simplify  the  symbols,  the  error  signals  are  still  denoted  as 

Ee,  E^,  Er  and  Ev,  respectively.  The  core  of  the  system  is  the 

optimum  state  estimater,  where  four  tracking  loops  intersect.  It 

could  complete  Kalman  filtering  or  other  forms  of  iterating 

filtering.  The  input  quantities  include  the*  four  error  signals 

Ec,  E^,  Er  and  Ev,  angular  velocity  components  u>e,  wd  and  wr  of 

the  electric  axis  of  the  antenna  revolving  with  respect  to  the 

inertia  coordinate  system,  antenna  control  commands  u>  and  «  . , 

ce  cd 

acceleration  components  ae,  a^  and  ar  of  the  target  in  the 
inertial  coordinate  system  as  given  by  the  inertia  guidance 
system,  and  two  directional  angles  and  0^  of  the  antenna 
coordinate  system  with  respect  to  the  carrier  coordinate  system. 

The  optimum  control  command  requires  the  optimum  state  estimater 

.  A  A 

to  output  estimated  values  of  angular  errors  Eg  and  E^  as  well  as 
angular  velocities  of  the  antenna  aiming  axis  w_  and  w_.  In  the 

E*  D  A 

meantime,  estimated  values  of  range  error  R  and  velocity  error  V 
are  also  given.  The  antenna  control  command  was  created  by  w,,=wi+Ae 
and  •  Simulation  results  showed  that  \  should  be 

A  /V 

10  ~  15.  In  addition,  R  and  V  were  used  to  form  range  gate  and 
velocity  gate  control  commands  R^c  and  VgC/  respectively. 

The  tracking  system  discussed  here  is  different  from 
classical  systems: 


Figure  1.  Structure  of  Tracking  System. 

Key:  1)  target;  2)  aiming  axis;  3)  electric  axis;  4)  antenna; 

5)  formation  of  range  gate  and  velocity  gate;  6)  formation  of 
range  and  velocity  commands;  7)  radar  receiver  and  signal 
processor;  8)  pre-balancer;  9)  optimum  state  estimater; 

10)  inertial  guidance  system;  11)  antenna  servo  mechanism; 

12)  antenna  command  formation. 

.  • 

(1)  A  Kalman  filter  was  used  to  complete  error  signal 

filtering,  i.e. ,  to  extract  the  true  estimated  value  of  an  error 
signal  from  noise  pollution  in  order  to  obtain  a  more  accurate 
deviation  signal  for  tracking. 

(2)  The  use  of  a  Kalman  filter  could  estimate  more  useful 
information  such  as  angular  velocities  «E  and  of  the  aiming 
axis  as  well  as  estimated  acceleration  components  a_.,  aL  and  a_ 
of  the  target  along  the  aiming  axis  coordinate  system.  These 
could  form  a  quasi-optimal  control  pattern. 

(3)  The  use  of  a  Kalman  filter  in  the  tracking  system  could 
naturally  form  a  hinge  for  all  tracking  loops.  A  tracking  loop 
could  extract  useful  information  from  other  loops  to  further 
improve  the  tracking  performance.  For  example,  it  was 
demonstrated  in  simulation  that,  when  the  range  and  velocity 
tracking  loops  intersected,  the  accuracy  in  ranging  could  be 
improved  from  the  classical  10  m  to  2  m  (mean  square  root) . 


(»)  R-i.  ■4800m  (bjRmi.  -800m 

Figure  2.  Comparison  of  Angular  Tracking  Error. 

Key:  1)  classical  tracking  system;  2)  this  tracking  system; 

3)  classical  tracking  system;  4)  this  tracking  system. 

(4)  The  parameters  in  the  filter  structure  are  time 
dependent.  Therefore,  the  optimal  tracking  performance  may  be 
maintained  over  a  very  wide  range  to  allow  the  system  to  have  a 
certain  adaptivity  to  target  maneuverability. 

In  order  to  compare  the  extent  of  improvement  of  the  Kalman 
filtering  tracking  system  relative  to  a  classical  system,  we  also 
conducted  the  simulation  of  a  classical  angular  tracking  system 
under  the  same  conditions.  The  sampling  frequency  of  the 
classical  simulation  system  was  103  Hz.  Results  of  the  /49 

simulation  are  shown  in  Figure  2. 


It  showed  that  the  tracking  system  discussed  in  this  paper 
had  good  steady  state  and  transient  characteristics.  The 
classical  system,  however,  required  a  longer  transition  time. 

The  steady  state  accuracy  could  still  be  maintained  for  far  away 
targets.  For  close  range  moving  targets,  the  tracking 
performance  deteriorated  apparently.  Thus,  the  tracking  system 


discussed  in  this  work  is  obviously  better  than  classical 
systems . 

II.  Target  Dynamics  Modeling 

Model  construction  is  the  key  to  applying  Kalman  filtering 
to  radar  tracking  systems.  For  an  airborne  tracking  system  for  a 
single  moving  target,  it  involves  the  inertial  coordinate  system 
(N,  E,  D) ,  carrier  coordinate  system  (X,  Y,  Z)  ,  antenna 
coordinate  system  (r,  e,  d)  and  aiming  axis  coordinate  system 

i 

(R,  E,  D) .  The  first  problem  is  which  coordinate  system  should 
be  used  for  target  dynamics  modeling. 

The  filtering  problem  in  a  radar  tracking  system  is  actually 
a  non-linear  problem.  In  engineering,  it  is  usually  linearized 
to  be  followed  by  an  estimation  of  Kalman  filtering.  It  is  a 
first  order  approximation.  Because  the  final*  solution  of 
non-linear  filtering  has  not  yet  been  obtained,  there  is  no 
rigorous  conclusion  on  the  convergence  of  approximations  of 
various  orders.  Therefore,,  we  evaluated  the  characteristics  of 
various  coordinate  system  models  through  simulation.  In  order  to 
obtain  linear  models  for  both  the  state  and  measurement,  we  also 
established  a  differential  polynomial  dynamic  model  in  the 
simulation  together  with  three  other  models  mentioned  earlier 
(carrier  coordinate  system,  antenna  coordinate  system,  and  aiming 
axis  coordinate  system) . 

The  same  typical  flight  pattern  was  used  for  the  four  models 
during  simulation.  In  the  first  4  seconds  the  relative  motion 
between  the  target  in  front  and  the  carrier  plane  was  linear  and 
turning  flight.  Later,  it  was  close  range  crossover  route 
flight.  Great  maneuverability  was  shown  in  this  case  [5] . 

Figure  3  is  a  series  of  simulated  curves.  Each  curve  was 
composed  by  110  data  points.  Results  showed  that  the  tracking 
accuracy  is  basically  the  same  for  tracking  systems  built  by 
these  four  models  when  the  mathematical  model  matches  with 
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Figure  3.  Comparison  of  Angular  Tracking  Accuracy  of  Four 
Models. 

Key:  1)  antenna  coordinate  system?  2)  aiming  axis  coordinate 

system;  3)  polynominal  model;  4)  carrier  coordinate  system. 

reality.  However,  under  maneuvering  conditions,  the  aiming  axis 
coordinate  system  and  differential  polynominal  dynamic  models  are 
more  adaptable  than  the  carrier  coordinate  system  and  antenna 
coordinate  system  models.  Therefore,  we  used  the  aiming  axis 
coordinate  system  in  this  work. 


The  aiming  axis  coordinate  system  was  used  in  the  /50 

airborne  radar.  The  first  problem  encountered  was  the 
transformation  from  an  inertial  coordinate  system  through  the 
carrier  coordinate  system  and  antenna  coordinate  system  to  the 
aiming  axis  coordinate  system.  Next,  we  should  also  consider 
that  these  four  coordinate  systems  are  always  in  relative  motion. 
It  will  involve  complex  dynamic  problems.  In  order  to  find* the 
equation  of  state  of  the  radar  tracking  system,  long  mathematical 


derivations  were  required.  This  work  could  be  foi  '  c 
scientific  research  reports  16-7]  and  only  the  fii.cx  r^si  re 

given  here.  The  equation  of  state  of  the  radar  tracking  syttem 
is: 
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Here,  the  first  order  Markov  process  is  used  to  describe  the 
target  acceleration,  i.e.,  &—  — (i/r)a+  wa  •  In  equation  (1),  r 
td/  tr*  ojd,  and  cjr  are  corresponding  components  of  r  and  w. 

III.  Linear  Filtering  and  Simulation 

Obviously,  equation  (1)  is  non-linear.  It  Cu,;  be  simply 
denoted  as 


J(D=/(*(0,0+G»(0  (2) 

It  is  difficult  to  realize  in  real  time  tracking.  In  order  to 
linearize  the  equation  of  state,  equation  (2)  was  expanded  into 
Taylor  series  in  the  neighborhood  of  x(t)  [6],  which  is  the 
estimated  value  of  state.  By  taking  the  first  order 
approximation,  the  linearized  equation  is  [6] 

i<*)=ACt)*(0+“C0  +  Cw(0  (3) 

The  difference  between  equations  (2)  and  (3)  lies  in  A(t)  and 
u(t).  Respectively,  A(t)  and  u(t)  in  equation  (3)  are  [6]: 
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It  should  be  noticed  that  A(t)  and  u(t)  are  functions  of  x(t)  and 
they  are  tine  dependent. 

t 

We  should  also  consider  the  effect  of  angular  blinking.  A 
first  order  Markov  process  could  be  used  to  describe  angular  /51 
blinking 

$<(0  =  “<1A*>s<0>  +  «',4(0 


Furthermore,  a  state  expansion  method  was  used  to  consider 
angular  blinking  as  two  state  variables.  The  effect  of  angular 
blinking  on  the  accuracy  of  the  estimation  could  then  be  reduced. 
Thus,  the  state  variables  would  be  expanded  to  11  dimensions,  not 
suitable  for  real  time  estimation.  Therefore,  we  considered  the 
angle,  range  and  velocity  channels  as  statistically  independent 
channels  for  separate  filtering  [6] : 
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Equations  (8) ,  (9)  and  (10)  are  continuous  time  equations  of 
state.  The  discrete  equation  of  state  and  equation  of 
measurement  are  as  follows  [7] : 

x(*  + 1)  =*(*+ 1,  *>*<%>  +  B(.h+  D«<*>  +  ( 1 1  > 

+  «'<*> 


where  v(k)  is  the  measuring  noise;  primarily  including  the 
thermal  noise  of  the  receiver,  clutter,  random  measuring  error, 
quantized  error  and  non-linear  effect  of  the  receiver  when  larger 
tracking  errors  exist. 

The  least  square  of  the  linearization  of  the  state  vector  of 
the  dynamic  system  was  estimated  by  Kalman  filtering; 


x(^+l|^)=^a+l.^jfW+5(^+1)u(^ 

«*+m+i)=K/t+ii*>+K<*+i>£z<*+l>-tf<*+1)K*+1|;c)]  .  (14) 

/ca+i)=pa+iU)//T^+i>c//^+i)P('c+H'c)«T('c+i>+i?^+I):i‘  <15> 

P(k  4-1  I  *••>=(*(■& +  l.fc)POc)$T(JC+l*Jc)+QW  (16) 


P(k+ 1 \k+  D=U-K(K+  D//Oc  +  l'llPOc+i  !*>[/-*(*+  »//<*+  D]T 

+K(.K+i'>R(k+DK'ca+l)  (17) 

In  order  to  investigate  the  performance  of  the  tracking 
system  after  Kalman  filtering,  a  series  of  simulations  was 
carried  out  on  a  FELIX  C-256  digital  computer.  The  principle  of 
the  simulated  tracking  system  is  shown  in  Figure  4.  During 
simulation,  the  computer  created  a  real  flight  condition.  It 
also  produced  engine  noise,  measuring  noise  and  random  noise  such 
as  angular  blinking.  In  addition,  it  created  observed  quantities 
by  the  receiver  and  signal  processor  and  formed  the  Kalman 
filter,  servo  control  loop  and  antenna  stabilization  loop.  One 
can  refer  to  Reference  [5]  regarding  the  stability  of  airborne 
antennas,  and  Reference  [7]  regarding  the  flow  chart  of  the 
computer  program. 

The  radar  carrier  and  the  target  were  in  relative  motion  in 

a  three-dimensional  space  in  simulation.  After  the  target  was 

captured,  the  carrier  aircraft  began  to  fly  horizontally  with  an 

2 

acceleration  of  30m/ s  from  a  constant  speed  horizontal  /52 

flight.  The  target  turned  45°  to  the  right  in  front  of  the 
carrier  and  began  to  climb  at  an  angular  speed  of  = 

\Zl50i+i50*  =212m  rad/s  •  In  addition,  it  exercised  evading 
maneuvers  along  the  positive  axes  of  the  aiming  axis  coordinate 
system  at  an  acceleration  of  2G.  Monte  Carlo  experiments  showed 
that  1/tr  =  1/TE  =  1/TD  =  1/20  during  the  evading  maneuver.  The 
step  used  in  Kalman  filtering  calculation  was  T  =  0.1s.  When  the 
control  command  was  created,  the  weighed  coefficient  X  =  15. 

Figure  5  shows  the  characteristic  curves  of  a  tracking 
system  using  Kalman  filtering.  EAch  curve  was  plotted  from  50 
data  points.  Here,  only  the  first  5  seconds  were  represented. 
Afterward,  it  was  completely  in  steady  state  tracking.  From 
simulation  one  can  see  that  the  radar  tracking  system  has 
relatively  higher  angular  and  range  tracking  accuracies  after 
adopting  Kalman  filtering.  Furthermore,  it  is  also  capable  of 


Figure  4.  Principle  of  Tracking  System  Simulation. 

Key:  1)  Relative  Motion  State  between  Target  and  Carrier; 

2)  angular  noise;  3)  formation  of  mixed  noise  from  errors  in 
angles,  range  and  velocity;  4)  range  and  velocity  noise; 

5)  optimal  filter;  6)  control  mechanism;  7)  antenna  stabilization 
loop;  8)  real  state  generator. 

accurately  estimating  the  angular  velocity.  •  It  is  adaptable  to 
tracking  maneuvering  targets. 

IV.  Linear  Filtering  Approximation 

The  aforementioned  linearization  treatment  should  be 
permissible  in  engineering  applications.  However,  the  gap 
between  this  approximation  and  the  optimum  should  be  known.  We 
should  also  understand  whether  there  is  a  need  to  reduce  this  gap 
and  the  expense  required  to  do  so.  To  this  end,  we  simulated 
non-linear  filtering  under  the  same  condition  as  linear  filtering 
and  obtained  the  appropriate  conclusions  through  comparison. 
Figure  6  shows  the  mean  square  root  error  curves  for  angle  and 
angular  velocity  tracking  [8] . 

It  indicated  that  the  use  of  non-linear  filtering  in  angle 
tracking  in  a  transient  process  was  slightly  better  than  Kalman 
filtering.  The  angle  error  was  improved  by  5-10%  and  the  angular 
velocity  error,  by  10-15%.  In  steady  state,  however,  there  was 
no  improvement.  We  also  observed  in  simulation  that  the  tracking 


Figure  5.  Performance  of  the  Tracking  System 

(a)  pitch  angle  tracking  error;  (b)  azimuthal  angle  tracking 
error;  (c)  pitch  angle  velocity  estimation;  (d)  azimuthal  angle 
velocity  estimation;  (e)  angular  tracking  mean  square  root  error 
(f)  range  tracking  mean  square  root  error. 


Figure  6 .  Mean  Square  Root  Tracking  Errors . 

(a)  mean  square  root  angle  error;  (b)  mean  square  root  angular 
velocity  error. 


Figure  7.  Effect  of  Noise  on  Tracking  Accuracy 

(a)  Effect  of  Measurement  Noise  on  ANgular  Tracking  (small 
dynamic  range) 

(b)  Effect  of  Noise  on  Tracking  Accuracy  of  Angles  (large  dynamic 
range) 

(c)  Effect  of  Noise  on  Estimation  of  Angular  Velocity  (large 
dynamic  range) 

Key:  1)'  model  in  agreement  with  reality;  2)  square  deviation  of 
engine  noise  increased  by  10  ;  3)  square  deviation  of  measurement 
noise  increased  by  5*;  4)  model  in  agreement  with  reality; 

5)  square  deviation  of  engine  noise  increased.by  10  ;  6)  square 
deviation  of  measurement  noise  increased  by  5  . 
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performance  in  range  and  velocity  was  not  improved  by  adopting 
non-linear  filtering  no  matter  whether  it  was  in  either  transient 
or  steady  state.  This  further  indicates  that  it  is  appropriate 
to  use  Kalman  filtering.  There  is  no  need  to  spend  several  times 
more  computations  in  non-linear  filtering  than  Kalman  filtering. 
Furthermore,  because  there  is  a  characteristic  "threshold”  in 
non-linear  filtering,  when  the  signal-to-noise  ratio  is  lower 
than  a  certain  threshold,  the  performance  may  drastically 
deteriorate.  It  is  not  favorable  in  all  conditions. 

V.  Sensitivity  Simulation 

The  optimal  condition  for  Kalman  filtering  is  that  the 
filter  model  coincides  with  the  actual  situation.  Generally,  it 
is  not  possible.  The  purpose  of  sensitivity  simulation  is  to 
analyze  the  effect  of  this  difference,  i.e.,  to  observe  the 
sensitivity  to  changes  in  engine  noise  and  measurement  noise  to 
determine  the  adaptability  of  the  filter  model. 

The  sensitivity  curves  of  the  tracking  system  discussed  /54 
in  this  work  were  obtained  through  simulation.  Figure  7(a) 
shows  the  angle  tracking  error  as  the  measurement  noise  varied  in 
a  small  dynamic  range.  Figure  7(b)  shows  the  effect  of  noise  on 
angle  tracking  in  a  large  dynamic  range.  Figure  7(c)  is  the 
effect  of  noise  on  the  estimated  angular  velocity  in  a  large 
dynamic  range. 

The  following  conclusions  are  obtained  from  the  simulation: 

(1)  When  the  square  deviation  of  the  measurement  noise 
varies  in  a  small  dynamic  range,  the  tracking  error  will  increase 
with  increasing  square  deviation.  However,  the  changes  are  not 
large  and  have  no  significant  effect  on  the  steady  state.  This 
system  is  not  sensitive  to  noise  changes  in  a  small  dynamic 
range . 


(2)  When  the  square  deviation  of  measurement  noise  varies  in 

2 

a  large  dynamic  range  (by  5  ) ,  the  tracking  accuracy  apparently 
decreases.  The  response  of  the  system  is  very  sensitive. 

(3)  When  the  square  deviation  of  engine  noise  varies  in  a 

2 

large  dynamic  range  (by  10  ) ,  the  decrease  in  the  tracking 
accuracy  is  far  less  than  that  due  to  measurement  noise,  i.e., 
not  as  sensitive. 

(4)  This  system  is  not  sensitive  to  engine  noise, 
demonstrating  its  maneuvering  tracking  capability.  This  is  what 
we  anticipate.  The  effect  of  measurement  noise,  however,  can  be 
corrected  by  measuring  the  signal  to  noise  ratio. 

VI.  Adaptive  Tracking  of  Maneuvering  Target 

When  Kalman  filtering  was  first  used,  maneuvering  was 
described  by  a  zero  average  white  Gaussian  process  which  is  easy 
to  treat  mathematically.  Obviously,  the  adaptivity  of  the  filter 
to  a  maneuvering  target  was  poor.  In  1970,  R.  A.  Singer  used  the 
zero  average  correlated  random  acceleration  to  describe  a  random 
maneuvering  [9] .  We  had  used  this  method  in  the  simulation 
mentioned  earlier.  However,  this  model  actually  assumes  that  the 
maneuvering  acceleration  is  a  steady  process.  Once  a  filter  is 
designed,  it  can  only  adapt  to  a  typical  maneuvering  pattern.  In 
order  to  overcome  the  shortcoming  of  the  aforementioned  method, 

R.  L.  Moose  proposed  a  model  in  which  a  target  maneuver  is 
decomposed  into  the  sum  of  white  noise  perturbation  and  finite 
maneuvering  commands  in  1973  and  1977.  In  1979,  R.  L.  Moose, 
et  al.  improved  the  model  to  a  composite  maneuver  model.  It  is 
even  closer  to  the  real  situation,  i.e.,  the  sum  of  the 
correlation  process  and  several  maneuvering  instructions  [10] . 

In  this  case,  the  adaptive  filter  was  comprised  of  a  Kalman 
filter  and  control  command  estimater.  The  semi-Markov  process 
was  used  to  describe  the  maneuvering  control  command.  For  the 


1. 


Figure  8.  Adaptive  Filtering 

(a)  Adaptive  Filter;  (b)  Angle  Tracking  Error  under  Stepwise 
Command 

Key:  1)  input  from  measurement;  2)  filter  I;  3)  filter  N; 

4)  filter  gain;  5)  gain  correction;  6)  weighed  probability; 

7)  adaptive;  8)  non-adaptive . 

first  time,  we  applied  this  model  to  the  aiming  axis  coordinate 
system  of  a  tracking  system  and  performed  a  simulation  [11] . 

The  adaptive  filter  system  is  a  series  of  parallel  filters 
as  shown  in  Figure  8(a).  These  filters  can  be  described  as  [11] 


i(t  lt> (*>£<(*  I  *>  (If 

i  •! 

Ultimately,  an  adaptive  filter  can  be  built  by  a  Kalman  / 5J 

filter  and  a  Bayes  estimator.  It  can  be  written  as  [11] : 


*«+ 1 \k+ o  -•<*+ 1.  +«tt) +m+ oc*«+ o 

-£ra+i>»a+i,«i«i*)+rtt)fitt)3  (19) 

where  u(k)  is  the  estimated  value  of  the  random  control  command: 

u(*)=ZX«'.<*+ 1)  (20) 

f  •  1 

The  weighted  probability  can  be  calculated  based  on  the  following 
formula  [11] : 

"  (21) 

where  W  (k) =Zwi(.k)  ,  c(k)  is  a  normalization  factor, 

and  A(k)  is  a  time-dependent  diagonal  matrix  whose  elements  can 
be  calculated  from  measurement  residue  and  input  command. 
i“  the  domain  transformation  probability  matrix  in  the  semi- 
Markov  process.  It  can  be  calculated  based  on  the  probability  of 
deviation  from  linearity.  The  covariance  matrix  of  the  estimated 
error  of  an  adaptive  filter  can  be  corrected  by  the  following 
formula  [11] : 

Pm  ik  I  K>  *  P(k  I  h) + {S  (kKii  ik  I K) -£(k  I  WCi,  (k  I  k)  -£(k  I  k)  JT  }  (22) 

Under  the  conditions  of  step  command  function.  Gauss  command 
function  and  exponent  command  function,  we  simulated  adaptive 
filtering.  Figure  8(b)  is  the  adaptive  and  non-adaptive  angle 
tracking  error  curves  under  a  step  random  command.  From 
simulations  one  can  see  that  the  composite  maneuvering  model  is 
superior  to  that  proposed  by  Singer.  In  a  transient  state,  the 
composite  maneuvering  model  is  25-35%  better  than  that  of 
Singer's.  But,  they  are  basically  identical  in  steady  state 
tracking. 


VII.  Conclusions 


(1)  The  Kalman  filtering  airborne  radar  tracking  system  is 
superior  to  classical  systems.  It  is  more  obvious  in  the 
presence  of  a  maneuvering  target. 

(2)  It  is  appropriate  to  use  the  aiming  axis  coordinate 
system  in  an  airborne  radar  tracking  system  for  a  single  target. 
As  the  coordinate  system  is  in  motion  it  facilitates  the 
linearization  of  the  model.  The  system  has  a  high  tracking 
accuracy  and  good  adaptivity  to  maneuverability. 

(3)  The  system  has  a  high  velocity  tracking  capability.  Th 
transient  tracking  characteristics  of  the  system  can  be  improved 
by  choosing  the  composite  maneuvering  model. 

(4)  Measurement  noise  is  one  of  the  major  factors  causing 
tracking  error.  In  a  real  system,  the  signal-to-noise  ratio  can 
be  introduced  in  the  calculation  to  determine  changes  in  the 
square  deviation  of  noise  in  order  to  correct  the  filtering  gain 

(5)  Tne  weighed  index  is  found  to  be  \  =  10~15  in  the 
simulation. 


Mao  Shiy ,  Li  Shaohong,  Zhan  Shuliang,  Wang  Guoying,  Zhang 
Xinghua  and  Zhang  Zhaowu  also  participated  in  this  work. 
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A  MODIFIED  KALMAN  FILTER  FOR  TRACKING 
MANEUVERING  TARGETS* 
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Fang  Ji-chang 

(Nanjing  Research  Institute  of  Electronic  Technology) 

Abstract:  A  modified  Kalman  filter  for  tracking  a 
maneuvering  target  was  presented  in  this  work  based  on  the  model 
of  filter  bank  described  in  reference  [1] .  It  detects  a  target 
maneuver  by  judging  whether  there  is  a  bias  in  the  observation 
residue.  The  maneuver  acceleration  command  is  only  estimated  by 
the  damping  least  square  method  when  target  maneuver  is  detected. 
Furthermore,  this  estimate  is  used  to  correct  the  state- 
prediction  and  error  covariance.  Otherwise,  it  works  as  a  single 
Kalman  filter  with  its  maneuver  acceleration  command  at  zero. 
Thus,  a  better  compromise- between  the  steady  state  filtering 
accuracy  and  fast  response  to  maneuver  is  reached.  Computer 
simulation  results  showed  that  the  accuracy  of  the  filter 
introduced  in  this  work  is  slightly  superior  to  that  of  the 
complicated  filter  bank  mentioned  in  Reference  [1] .  The 
computation  load,  however,  was  only  1/3.6  of  the  latter  model. 

Abstract:  A  Modified  Kalman  filter, bated  on  the  model  of  filter  bank  for  tracking  ma¬ 
neuvering  targettftee  referenced 31  it  pretented  inthitpeper.lt  detect!  the  target  maneuver  bp 
judging  whether  obtervation  retidue  exhibit!  biat.  The  maneuver  acceleration  command  is  es¬ 
timated  bp  damp  least  squares  method  only  when  the  maneuver  it  detected,  .and^^n  the  filter 
will  correct  the  state-prediction  and  its  error  covariance  with  the  estimate.  Otherwise  the  fil¬ 
ter  works  at  a  single  Kalman  filter  with  zero  maneuver  acceleration  command.  Thus,a  better 
compromise  between  the  steady  filtering  accuracy  and  rapid  response  to  maneuver  is  obtained. 

The  computer  simulation  results  show  that  the  filtering  accuracy  is  somewhat  better  than  that 
of  Moose's  complicated  filter  bank  and  the  computation  burden  is  only  1/3.6  of  Moose’s  filter. 

I .  Introduction 

In  Reference  [1] ,  the  acceleration  motion  of  a  maneuvering 
target  was  considered  as  a  combination  of  a  semi-Markov  process 
[2,  3]  and  the  Singer  model  [4] .  A  filter  bank  comprosed  of  N 
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Kalman  filters  was  used  to  find  N  possible  filtering  states  and 
their  corresponding  probabilities  with  respect  to  N  possible 
maneuver  acceleration  commands  u^,  ieN.  Then,  their  weighted 
average  was  obtained  as  the  state  filtering  value  at  that  precise 
moment.  The  mathematical  model  of  this  method  is  more  accurate. 
Under  the  premise  that  the  target  maneuvering  characteristics  are 
known  in  advance,  the  filtering  accuracy  is  relatively  higher. 
However,  the  computation  load  is  also  heavy.  It  is  also  highly 
dependent  upon  prior  knowledge.  Therefore,  the  assumptions  made 
for  may  not  be  accurate  in  actual  applications.  In  addition, 
all  the  main  diagonal  elements  in  the  Mafkov  probability 
transformation  matrix  are  P  and  the  remaining  ones  are  (l-P)/(N-l) 
and  (P—1).  This  assumption  also  does  not  completely  agree  with 
the  actual  situation.  Therefore,  the  actual  filtering  accuracy 
decreases. 

Most  of  the  targets  considered  in  this  work  were  non¬ 
maneuvering  within  relatively  long  periods  of  time.  They  were 
maneuvering  only  small  portions  of  time.  They  exhibited  strong 
maneuvers  only  in  even  smaller  fractions  of  time.  Hence  the 
modified  Kalman  filter  was  comprised  of  a  single  u^  =  0.  Singer 
type  Kalman  filter  with  a  maneuver  detector  and  a  damp  type  least 
square  estimater  for  maneuver  acceleration  command.  /58 

Consequently,  in  a  relatively  long  period  of  time  in  which  u^  is 
indeed  zero,  the  model  coincides  with  the  reality.  The 
computation  load,  however,  is  lower  than  that  for  a  filter  bank. 
Furthermore,  because  the  observation  residue  is  used  to  correct 
the  model  while  filtering,  the  deviation  caused  by  u^  M  can 
also  be  better  rectified. 

II.  Model  Analysis 

The  equation  of  state  and  equation  of  measurement  for  the 
system  given  in  Reference  [1]  are: 


(1) 


z,(^+i)=H,a+i)x,a+i)+v,a+i) 

where  j  =  1,  2,  3,  representing  the  three  channels-range  r,  angle 
of  elevation  e  and  orientation  b,  respectively.  VT  (k)  and  (k) 
represents  model  noise  and  measurement  noise,  respectively.  They 
are  both  zero  average  independent  white  Gauss  noises,  u^ (k)  is 
the  maneuver  acceleration  command.  Let  us  use  the  range  channel 
as  an  example,  Xi(t+1),  A),  (AiC^  +  lA)»  Wi(^)  are 

expressed  by  the  following  formula:  7 
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whereAsd-tf—O/o,  fl=<l  +  <a«-"-aa-")/<a-a)]/(aa),  C-  (oT- 1  +  «-•*)/ 
a*,  £>=Cr+(aA-o/)/(o-d):/(aa)«,  E=e‘aT,  F  =  (ff-*T-£)/(a -a)  ,  G  =  (/~ 
A)/(a-a),  /=(1 /a,  a=0.4/s,  a=0.l/s. 


(3) 


T  is  the  sampling  time  interval.  Z1 (k  +  1)  is  the  measured  value 
of  range  r(k  +  1).  When  only  the  range  is  measured,  H^(k  +  1)  = 
[100]  . 


Let  =  0 .  The  single  channel  Kalman  filtering  equations, 
i.e.,  the  expressions  for  state  prediction,  state  filtering 
value,  gain,  covariance  of  state  prediction  error  and  covariance 
of  state  filtering  error  are: 


=^(-^-)+KU+l)[z(/c+l)  -«0c+l>x(-^p-)] 

/ca+i)=p(-^-)HTa+i)  [h(.k+\'>p(-^jl-)ht(.k+i)+r(,k+i')]  ‘ 

p(l$r)=  CI-KOt+DHa+DJpf-^j1-) 


(4) 


*In  Reference  [1],  o-CT  +  <aA-aT)/«i-a)]<aa),  .  it  is  believed  to  be  mistaken. 
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where  Q(k)  =  E[W^(k)],  R ( k)  =  E[V^(k)],  and  I  is  a  unit  matrix. 

When  u^  #  0  ,  the  following  problems  will  arise  if  equation 
(4)  is  still  used  for  filtering.  1.  The  state  prediction  error 
X(k  +  1/k)  and  observation  residue  Z(k  +  1/k)  become  larger. 
Furthermore,  both  are  biased,  because 

1+  (-^r-)3?  (t)+#'  (JTJ-),v'«>+r(-TJ-)“«)> 


ua>=o 
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Correspondingly,  L  '  | ^ u(.K)^0 

Then,  2a+iA)=za+l)~HC/c+i)^a+lA)=H(/c+l)^a+iA)+V(/c+i),  ( 

fO,  u(^)=0 

e\  zC^r-)}  =H(k+l)E[%  =  a+M  <56) 

L  v  K  /J  i  \  h  yj  o 


Here,  X(k/k)  =  X(k)  -  X(k/k),  which  is  the  state  filtering  error. 
Because  it  is  assumed  that  u(k)  #  0  beginning  at  kT,  therefore, 
E(X(k/k)]  =  0  is  still  valid.  2.  The  covariance  of  state 
prediction  error  becomes  correspondingly  larger: 
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where  Qu(k)  =  u  (k)  .  The  lower  equation  in  (5c)  is  an 
approximation.  One  can  see  that  it  is  necessary  to  correct  the 
state  prediction  and  error  covariance  of  a  filter  assuming  u(k) 
u^  =  0  when  u(k)  ^  0.  This  will  require  us  to  measure  and 
estimate  u(k). 

III.  Maneuvering  Target  Detector 

In  Reference  [1]  the  observation  residue  Z(k  +  1/k)  was 
approximated  as  a  Gaussian  process.  Its  jsquare  deviation  is 

p.tt+i/«=£Cza+iA)zTa+iA)]=Va+i)pa+iA)HTa+i)+Ra+i). 

and  its  mean  value  is  EZ2ik+  1A)1  •  In  fact  equation  (5b)  is 

different  when  u(k)  =  0  and  £  0.  This  is  the  basis  for  using  Z 
as  a  statistical  check  for  u(k) .  As  shown  in  Figure  1,  |zj in  the 
figure  should  be  Z. 

For  this  reason,  a  coincidence  detector*  should  be  used.  It 
criterion  is 

I21>*VpT  (6) 

If  it  is  satisfied  consecutively  n  times,  then  it  is  determined 
that  there  is  a  maneuver  command.  Otherwise,  it  is  determined 
that  a  maneuver  command  does  not  exist.  Obviously,  when  m  and  n 
increase,  the  false  alarm  probability  and  spotting  probability 
also  decrease  correspondingly.  Otherwise,  they  increase 
accordingly.  Here,  «=2. 


Figure  1.  Probability  Density  Curve  of  Observation  Residue. 


IV.  Damped  Least  Square  Estimation 


When  u(k)  ^  0,  it  is  also  necessary  to  calculate  u(k)  to 
make  the  correction.  The  relationship  between  u(k)  and  the 
observation  residue  must  be  derived  first  and  then  u(k)  can  be 
obtained  using  the  damped  least  square  method. 

When  u(k)  #  0,  the  state  estimation  based  on  filters  where 
u(k)  =  0  is  biased  (expressed  with  a  superscript  "  '  ").  The 
symbol  for  the  unbiased  estimation  remains  unchanged.  The  symbol 
for  the  corresponding  observation  residue  also  remains  the  same. 
It  is  not  difficult  to  prove  that  the  prediction  of  biased 
estimation  and  its  observation  residue  are:  /60 
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The  unbiased  estimation  and  its  observation  residue  are: 


2('^)=Za+1)‘Ha+1)^(jL^L)  (8 

z(-|$f)=z«+2>“««+2)^(^r) 

After  subtracting  equation  (7)  from  equation  (8)  and  moving  the 
terms,  one  can  get: 
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If  the  duration  of  u,  t,  satisfies  r»l7\  then  u(*)«u(^+l)ss  /  61 
•••ssuOt+O^u'  •  Hence,  equation  (9)  may  be  written  in  a 

matrix  form 


Z'  «Cu'  +Z 


(10) 


One  can  obtain  the  following  by  using  the  weighted  least  square 
estimation  (Markov  estimation)  [5] : 


u=[CTp;1cr,CTP;1z/,  p.=ezzzti.  (ID 

The  more  data  points  there  are  in  equation  (9) ,  the  higher  the 
accuracy  becomes.  However,  the  computation  load  and  memory 
requirement  are  also  higher.  Furthermore,  there  is  more  delay 
involved.  In  order  to  synchronize  real  time  processing  and  the 
aforementioned  maneuver  detection  process,  two  points  were  used 
in  a  trial  computation,  i.e.,  to  take  £'  and  C  and  the  2 
dimensions  involved.  We  get 


,  CuPh(K+2/k+l)2'(k+l/k)  +  CnP.<k+i/k)l!,Oi+2/k+i'i 
u  ~  ciPAK+2/h+iy+ciP.cK+i/h) 


(12) 


When  equation  (10)  was  used  to  estimate  u',  the  component  of 

became  smaller  as  T  decreased.  The  situation  that  the  condition 
T 

C  C  deteriorates  as  the  least  square  method  is  used  to  solve  the 

contradictory  equation  begins  to  appear.  Or,  Cu'/ 

is  considerably  smaller  than  the  signal  to  noise  ratio,  affecting 

the  accuracy  of  estimation.  In  this  case,  the  damped  least 

square  method  may  be  used.  This  method  involves  adding  an 

appropriately  small  positive  number  Y  to  the  major  diagonal  line 
T 

of  C  C.  Y  is  continuously  iterated  and  varied  until  the 
satisfactory  accuracy  is  obtained  [6] .  It  is  not  suitable  to  use 
an  iterative  method  in  real  time  analysis’.  Therefore,  we  let 
CTP;1C+y/=pCTP;1C,  .  Furthermore,  the  optional  p  determined 
experimentally  ahead  of  time  is  used  in  equation  (11)  .  We  get 

5=CPCrP;1a-1CTP;1Z/=u,/p  (13) 

p  varies  with  T  (as  shown  in  Table  1) .  Its  empirical  formula  is 
as  the  following: 


P.=Cr,/ru,  rasCoT-l+e-^/a’j 

P»  =  1.2P„  P,  =  Cr,(l-e-“T)P.  J  (14> 


When  u(k)  #  0 ,  u  is  detected  and  estimated  to  perform  the 
following  correction:  ,  .  ,  .  „ 


(15) 


(15) 


+r  (±V~)i'r'  ("^TT") 


(16) 


Nc> 


Table  1.  Coefficients  CTg  and  CTr  Corresponding  to  Sampling 
Interval  T  and  t>e,  ofa. 


T 

Crr 

fir 

fir 

fit 

0.5*. 

1.17 

1.93 

3.5 

10 

12 

It 

1.27 

1.92 

1.9 

2.9 

3.5 

V.  Computer  Simulation  Results 

Parameters  in  Reference  [1]  were  used  to  generate  the  target 
trajectory  data  U|s{±l4g,  ±7g,  0>, .  Approximately  35%  of  the  time 
was  used  as  the  target.  The  mean  square  root  of  model  error  is 
J*=<T»s=flr*“l5n,/sI  The  mean  square  root  of  measurement  error  is 
<r,=25m,  a,  =at=3mrad.  Both  errors  were  generated  by  a  pseudo-random 
number  program.  The  measurement  error  was  larger  than  that 
assumed  in  reference  [1]  in  order  to  examine  the  performance  of 
the  filter  at  a  lower  radar  measurement  accuracy. 

The  target  tracking  data  generated  was  processed  by 
filtering  methods  described  in  the  paper  and  Reference  [1]  , 
respectively.  The  filtering  values  obtained  such  as  r>  f.  <•„  j  y 
and  their  corresponding  computer  processed  data  and  symbols  were 
plotted  as  shown  in  Figures  2(a)  -v  (d)  (only  a  series  of 
examples) .  Table  2  also  compared  the  filtering  value  of  this 
filter  and  the  Moose  filter  against  the  actual  number  with 
respect  to  T  *  0.5,  0.75,  1,  and  1.5s.  A  is  the  maximum 

C 

filtering  error  of  this  filter  and  A^  is  the  maximum  filtering 
error  of  the  Moose  filter  bank.  N  represents  the  number  of  times 
the  value  of  this  filter  comes  closer  to  the  actual  number  than 
that  of  the  Moose  filter  out  of  100  points  after  cancelling  the 
number  of  times  the  two  filter  values  approach  the  true  value. 
From  Table  2  one  can  see  that  in  most  comparisons  N  is  a  positive 
number.  This  means  that  there  are  more  times  this  filter  came 
closer  to  the  actual  value  than  those  of  the  Moose  filter.  / 6 2 

Its  maximum  error  is  also  slightly  smaller  than  that  of  the  Moose 


3. 

I#*  r-i.s 


r  »o.rs* 


Figure  2.  Results  of  Simulations  of  Two  Filters. 

-  represents  real  value 

*  represents  filtering  values  of  this  filter 

+  represents  filtering  values  of  the  Moose  filter 

.  • 

Key:  1)  range;  2)  orientation;  3)  angle  of  elevation;  4)  range 
speed. 

filter.  Hence,  the  accuracy  of  this  filter  is  slightly  higher 
than  that  of  the  Moose  filter  bank. 


-  w  •<*•:•**  >  1 


In  computer  simulation,  the  time  required  for  each  filtering 
method  was  gathered  and  printed.  On  the  same  computer,  the 
filtering  of  the  same  100  data  points  requires  57.96s  using  the 
Moose  filter.  It  only  takes  16.02s  using  this  filter.  The 
latter  is  1/3.6  of  the  former.  Obviously,  the  computation  load 
of  the  filter  introduced  in  this  work  is  significantly  reduced. 
Therefore,  there  are  reasons  to  believe  that  this  filter  is  more 
suitable  for  real  time  processing  than  the  Moose  filter. 


v\-'  Vi 


VI.  Conclusions 


The  maneuvering  target  tracking  filter  presented  in  this 
paper  uses  observation  residue  as  a  check  to  measure  and  estimate 
the  unknown  model  parameter  -  maneuver  acceleration  command.  The 
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Table  2.  Comparison  Chart  of  Accuracies  of  Two  Filters 
Key:  1)  compared  quantity. 
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function  of  the  filter  matches  with  the  model  parameters  so  that 
the  computing  capability  is  reasonably  utilized.  Consequently, 
from  a  certain  sense,  it  is  an  "identifying  and  correcting  while 
filtering  model"  modified  Kalman  filter. 


» 


This  work  was  completed  under  the  guidance  of  Senior 
Engineer  Zhang  Zhizong  and  comrade  Qian  Huisheng.  The  author 
wishes  to  express  his  gratitude. 


Vol.  AES-15,  pp.  448—456, 


C  I  3  R.  L.  Mom*,  H.  F.  Vanlandin|han  and  D.  H.  Meeabt,  IEEE  Tran* 

Map.  107V. 

CJ3  N.  H.  OhoJaoa  and  R.  L.  Moo...  IEEE  Tran..,  Vol.  AES-13,  pp.  310-317,  Map  1077. 
1.33  R.  A.  Howard,  IEEE  Tran*.,  Vol.  MIL-8,  pp.  114—124,  A| 1064 
C43  *•  A.  Sin(ar,  IEEE  Trana.  Vol.  AES-6,  pp.  473-483,  July  ift70. 

£53  axe.  aae-H-ac*.  ±*2a**,  1070*6,8. 

c«3  mmtxm,  «*«««,  107a*. 


vy 


END 


FILMED 


1-85 


,  ■  .  j 

5>:i 


DTIC 


•  r  •.•  .  •.- 


